Electric field and exciton structure in CdSe nanocrystals 
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Quantum Stark effect in semiconductor nanocrystals is theoretically investigated, using the ef- 
fective mass formalism within a 4 x 4 Baldereschi-Lipari Hamiltonian model for the hole states. 
General expressions are reported for the hole eigenfunctions at zero electric field. Electron and hole 
single particle energies as functions of the electric field (Eq_d) are reported. Stark shift and bind- 
ing energy of the excitonic levels are obtained by full diagonalization of the correlated electron-hole 
Hamiltonian in presence of the external field. Particularly, the structure of the lower excitonic states 
and their symmetry properties in CdSe nanocrystals are studied. It is found that the dependence of 
the exciton binding energy upon the applied field is strongly reduced for small quantum dot radius. 
Optical selection rules for absorption and luminescence are obtained. The electric-field induced 
quenching of the optical spectra as a function of Eqd is studied in terms of the exciton dipole 
matrix element. It is predicted that photoluminescence spectra present anomalous field dependence 
of the emission lines. These results agree in magnitude with experimental observation and with the 
main features of photoluminescence experiments in nanostructures. 

PACS numbers: 73.21. La,73.22.-f, 78. 67.-n 
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I. INTRODUCTION 



Semiconductor nanostructures under longitudinal 
electric field produce pronounced effects on optical 
propertiesiiiSiMiSiSii It has been shown that the field 
induces a red shift of the exciton peaks in the photo- 
luminescence and electro-optical spectra. The shift of 
the excitonic peaks to lower energy with the increasing 
electric fields is known as quantum Stark effect, while 
the reduction of the overlapping between the electron- 
hole pair wave function by the field is related to quench- 
ing of the fundamental transition in the luminescence 
spectrum. Zero dimensional systems as colloidal semi- 
conductor quantum dots (QDs) under electric fields are 
appropriate candidates for several device applications, 
including optical computing and fiber-optical commu- 
nication (see Ref. and references therein). Also, the 
micro-photolumincsccncc spectroscopic technique in sin- 
gle spherical QDs has allowed to study fundamental is- 
sues of the excitonic states^ 

The electro-optical properties, the Stark shift, and the 
dependence of the peak intensity in the optical spectra 
upon the applied field should depend strongly upon the 
details of the band structure. This was demonstrated 
for quantum wells in the 80's>i2iAiii2 It is interesting to 
investigate the analogous effects in QD, as the three- 
dimensional confinement causes properties that are be- 
yond the naive enhancement of the effects observed in 
quantum wells. For example, in QDs the band dispersion 
no longer exists, the energy spectrum is totally discrete 
and depends qualitatively upon the dot geometry. More- 
over, the high surface to volume ratio originates effects 
that are intrinsic to QDs. Perhaps, the most spectacular 
finding up to date is the discovery of high luminescence 
in porous Si^ where QDs are believed to play an im- 
portant role>i^ Other striking effects can be found in the 



dark magneto-exciton luminescence^ and the blinking 
and spectral shifting of single QD luminescence under 
external electric fields £ 

Calculations of the quantum Stark effect in spherical 
QDs in the strong confinement regime have been per- 
formed in the framework of the parabolic modelii£*i£ The 
simple parabolic model was able to provide a relative 
good picture for the description of the electronic states 
in the conduction band. This approximation breaks down 
for the calculation of the hole levels due to the fourfold 
degeneracy and the admixture of the light and heavy- 
hole bands present in the II- VI and II-III compounds 
with zinc-blende lattice structure. A reliable description 
of the energy band dispersion is ruled by the Baldereschi- 
Lipari Hamiltonian^ Within this approach, calculations 
of the influence of an external electric field Eqd were 
done in Ref. iPJ However, this calculation presents sev- 
eral limitations (for a detailed discussion see Ref. l20f) . As 
it is well known, the dependence of the interband optical 
transitions upon the light frequency (absorbed or emit- 
ted) reflects the structure of the conduction and valence 
bands. 

In this paper we study the excitonic Stark effect of 
spherical QDs taking into account the valence band 
admixture using the Baldereschi-Lipari Hamiltonian 18 . 
Supported by a rigorous treatment of the exciton wave- 
functions, we have obtained the interband dipole matrix 
elements taking into account the fundamental symme- 
try properties of the QD. From the dipole matrix ele- 
ments we have obtained the optical selection rules that 
allow the identification of the exciton levels observed in 
absorption and luminescence experiments. We present 
numerical calculations that reveal the combined effects 
of band admixture, Coulomb interaction, QD size and 
electric field intensity. 

In Sec. II we examine the energy dependence upon 
Eqb for the electrons and holes in CdSe QDs. Explicit 
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analytical solutions for the hole levels at Eq£> = are de- 
rived. The influence of Coulomb correlation and valence 
band coupling on the quantum Stark effect is analyzed in 
Sec. III. Section IV is devoted to study the electric-field 
induced optical properties. The main results of the paper 
are summarized in Sec. V. 



II. SINGLE PARTICLE STATES 

In the effective mass approximation, the Stark effect 
on the electronic states at the bottom of the conduction 
band can be represented by products of the T 6 Bloch 
functions (v\l/2,s z ) (with s being the conduction band- 
edge angular momentum and s z = ±1/2) times enve- 
lope functions. The later ones are obtained from the 
effective mass Hamiltonian with a uniform electric field 
\s\EQjjr cos 9. Here, e is the electron charge e and Eqd 
is the electric field intensity inside the nanocrystal. At 
zero electric field, the envelope functions take the form 
Rni{r)Yu 2 (9, (f). R n i(r) are the radial wave functions^ 
and Yu z (9, ip) are the spherical harmonics^ These states 
are in the I — s coupling scheme and have well defined 
values of the squared orbital and spin angular momenta. 
Instead, we will use the / — / coupling scheme, where 
the states have well defined total (f = 1 + s) angular mo- 
mentum projection hf z and square value h 2 /(/ + 1), that 
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where {l^l z s z \f fz) a re the Clebsch-Gordan coefficients, 
f = l z + s z , and \l - 1/2| < f < 1 + 1/2.21 If the hole 
states in the valence bands are described by an spherical 
4 x 4 k • p Hamiltonian, the I — / coupling scheme for 
the electron wave function is more convenient to build 
up the excitonic states in a spherical QD. The hole and 
electron states present the same symmetry properties, 
allowing to use all inherent properties of the Clebsch- 
Gordan coefficients. 24 

Since the s z = ±1/2 bands are uncoupled, the states 
described by Eq. Q have energies E^ l that are indepen- 
dent of the quantum numbers /, f z . In the simulation 
of the real electronic states, the confinement potential is 
chosen as a spherical box with an effective radius R e f, 
which is greater than the structural nanocrystal radius 
R. This effective radius is introduced in order to take 
into account, approximately, the penetration of the elec- 
tron wave function in the surrounding medium. In our 
calculations, we determine R e f from the condition that 
the energy of the Is state Ei$ = h 2 ir 2 /2m e R 2 ^ be equal 
to the energy calculated for a spherical well with depth 
V e = 600 meV& 

The electron states under external electric fields are 
found by numerical diagonalization of H e in the basis 
provided by The matrix elements of the Stark term 
are provided in the Appendix. 



TABLE I: Parameters used in the calculations. 
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For the hole states in the bottom of the va- 
lence band we use the well known Baldereschi-Lipari 
Hamiltoniani2*22i2£ in presence of a constant electric 
field, that is 



H h = 



7i 



2mo 
+V h (r) 



p2_g(p(2).j(2; 
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QD 



r cos ( 
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Vh(r) being the confinement potential for the valence 
band, and J< 2 ) are spherical tensors of rank 2 

built from linear and angular momentum operators, fi = 
272/71, and 72 and 71 are the Luttinger parameters of 
CdSe in the spherical approximation 72 = 73. The hole 
eigenfunctions at zero electric field can be cast as 
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£ £ (KIL Z J Z \FF Z ) 

K=L,L+2 L g ,J z 

xR { p K {r)Y KLz {9^) (r|3/2,J z ), (3) 



where (r|3/2, J z ) are the hole Bloch functions of the Tg 
valence band with band-edge angular momentum J = 
3/2. The hole Bloch functions are related with the va- 
lence electron Bloch functions \J,J Z ) by the rule \ JJ Z ) = 
{—l) J ~ Jz J J, — J z ) (derived from the time-reversal oper- 
ation). Our Bloch functions 1 3/2, J z ) have the following 
phase convention: |3/2,±3/2) = ^(i/y/2)(X ± iY) |±), 
and |3/2, ±1/2) = (i/%/6) [2Z |±) =F (X ± iY) |=f)]. The 
phase factors of the above Bloch functions arc implicit in 
the optical dipole matrix elements. 

For F = 1/2, according to the rule of the addition of 
two angular momenta, \L — 3/2 1 < F < L + 3/2, the 
states defined in Eq. 10 reduce to two uncoupled states 
L = 1 (K = 1) and L = 2 (K = 2) which correspond to 

(1/2) 
N,L 

ful- 



P1/2 and D 1 / 2 states with radial wave functions R 
[L = 1,2). In this case the eigenfunctions R 



(1/2) 
N,L 



fill two independent radial effective mass equations with 
light hole effective mass mih = mo/ (71 + 272). 

(F) 

For F > 3/2, the radial wave functions R N K (r) are 
solutions of the coupled differential equations 29,30 
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where W(r, E) = 2m [V h (r) - E] /h 2 ^. 

The coefficients Ou C 2 and C 3 are reported for sev- 
eral states in Refs. and |2|1 Following the argument 
of Baldereschi and Lipari, we have obtained the general 
expressions 



C X {L,F) = nVZ(-l) 



3/2+L+F 



L L 2 

3/2 3/2 F 



/2L(2i + l)(2L + 2) 



C 2 (L,F) = M V30(-l) 3/2+L+F 



(2L + 3)(2L-1) ' 

L + 2 L 2 
3/2 3/2 F 
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(L + l)(L + 2) 



C 3 (L.F) 



2L + 3 

d(L + 2,F) = -Cx(L,F). 



(5b) 
(5c) 



The result C3 = — C\ has been verified numerically. We 
also found numerically that Cf + C 2 = /s 2 and C 2 / V > 0. 

In the case of abrupt infinite confinement potential, 
the solutions of equations (0J are 



R*k> L (r)=A[j L (kr/R)- 



R 



{F) (r) 

N.L+2V I 

x (j L+2 (kr/R) 



A (c x /C 2 + V(Ci/C 2 ) 2 + 1 
j L +2(VPkr/R)j L+2 (k) 



,(6a) 



(6b) 



where (3 = mi h /m hh (m hh = m /(7i - 272) is the heavy 
hole mass. The parameter k fulfills the transcendental 
equation, 



3L+2 



(k) (<VC 2 



1 



2 1- V(Ci/c 2 ) 2 

(V0k)j L (k) (Cx/C 2 - V(Ci/C 2 ) 2 + l) , (7) 



and A is a normalization constant, such that 



R 



(F) 

AT, L+2 



l dr = 1. 



(8) 



The hole energies are equal to Eh = h 2 k 2 /2nihhR 2 ■ 

As in the case of the conduction band levels, the hole 
states under external electric field are found by numer- 
ical diagonalization of the Hamiltonian Hh in the basis 
provided by Eq. ©. The matrix element of the Stark 
term are provided in the Appendix. 

Figure n shows the energy levels of electrons and holes 
as functions of the electric field intensity. The energies 



are in units of Eq = ti 2 /2miR 2 {i = e, hh), and the elec- 
tric field intensity is in units of Eo/\e\R. In this figure, R 
for the electrons is the effective radius above introduced. 
The states at zero field are indicated by the usual spec- 
troscopy notation NA, with A = S, P, D, ...(s,p, d, ...) for 
the L = 0, 1, 2, ...(I = 0, 1, 2...) hole (electron) states. In 
the case of the hole energies, the quantum number F is 
indicated by a sub- index, that is, NAp- In the figure the 
values of \l z \ and \F Z \ for each states are indicated and, 
as it can be seen, at non-zero electric field the ± \l z \ and 
the ± \F Z \ degeneracy remains. For the electrons, due to 
the absence of Bloch-envelope (spin-orbit-like) coupling, 
there is an additional degeneracy. Thus, the electron en- 
ergy levels only depend upon the modulus of the angular 
momentum projection being degenerate in the spin 
projection. The results for the electrons here shown re- 
produce those of Ref. 0. The light hole energies, which 
correspond to F = 1/2, are higher and are not included 
in Fig. 1. 

The hole Hamiltonian (J2J does not include band warp- 
ing terms that arise from the cubic symmetry of the 
nanocrystal lattice. The cubic corrections are propor- 
tional to the cubic coupling parameter^ S = (73 — 72 ) /71 ■ 
The Luttinger parameter 73 has not been determined 
for CdSe, but it has been estimated as 0.532^ assuming 
that the ratio 73/72 is equal in CdTe and CdSe. Hence, 
a value S = 0.066 is obtained, which is one order of 
magnitude smaller than the spherical coupling param- 
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FIG. 1: Single particle energies vs. electric field intensity. 
Note that Eq = h 2 /2rriiR 2 depends upon the quasiparticle 
effective mass. 
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eter jj,. Therefore, the cubic corrections can be included 
using perturbation theory. For zinc-blende semiconduc- 
tors, parity-breaking terms exist, in principle. However, 
this effect seems to be smaller and the optical properties 
of acceptor levels have been explained without consider- 
ing them^ 1 For zero electric field, the hole states P\i% 
and -D1/2 belong to the irreducible representations Tq 
and F7 of the point group T^, respectively, while S3/2 
and P3/2 states belong to the irreducible representation 
r 8 . Hence, their energies do not split and the selection 
rules are not modified for these states. However, the 
energy values shift in second order of 5, unless coupled 
quasidegenerate states are present^ This is the case of 
the coupled states 26*3/2 — I.D5/2, where the splitting is 
first order in S. For each F > 3/2, the eigenstates of 
the cubic Hamiltonian are formed of linear combinations 
of states with different F z . These linear combinations 
transform according to the irreducible representations of 
T4 that are compatible with the representation Dp of the 
3-dimensional rotation group 0(3). Hence, the analyzed 
level splits proportionally to S in several levels, according 
to the compatibility tables of groups Td and <3(3). 32 The 
optical selection rules are relaxed for these states and 
additional transitions should appear with low intensity. 
The effects of the cubic anisotropy has been studied nu- 
merically for CdSe QDs using the Effective Bond Orbital 
Method (EBOM), 25 which implicitly takes into account 
the lack of inversion symmetry. It was shown that the 
spherical approximation is very good, specially for large 
nanocrystals. 

III. EXCITON STATES 
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FIG. 2: Energy levels vs electric field of the lower excitonic 
states of a CdSe nanocrystal 2 nm in radius. The inset shows 
the allowed optical dipole matrix elements squared (see Sec. 
II VI . By symmetry the components \D y \ 2 and \D X \ 2 are equal. 

The exciton Hamiltonian can be written as the sum 



of the electron and hole Hamiltonians plus the screened 
Coulomb interaction V e -h — — e 2 /e |r e — r/, |, e being 
the dielectric constant. Excitonic states, can be ob- 
tained using an expansion in a basis of electron-hole pair 
wave functions with well defined total angular momen- 
tum square h 2 M(M + I) and projection KM Z 

|*) = J2 C a \nNlLfF;MM z ), (9) 

where 

\nNlLfF; MM Z ) 
= UFf z F z \MM z )\nlff z )®\NLFF z ). (10) 

It is important to remark that in Eqs. © and (|10fl it 
is implicit the condition for the addition of two angular 
momenta for electrons and holes \ f — l/2|<Z</ + l/2 
and \F - 3/2| < L < F + 3/2, respectively. The matrix 
elements of the Coulomb interaction in the basis I© are 
reported in Ref. |33 The matrix elements of the Stark 
term are provided in the Appendix. 

As the electric field is chosen along the Z-axis, the exci- 
ton angular momentum Z-projection is a constant of mo- 
tion and the Hamiltonian can be diagonalized indepen- 
dently in different M z -subspaces. We have built the M z - 
subspaces using as basis all the possible electron-holes 
states HOI) that fulfill the condition 2m l R 2 E/h 2 < 90 
(i = e, hh). With these criteria the dimensions of the di- 
agonalized matrices are 532, 502, and 415 for \M Z \ = 0, 
1, and 2, respectively. 

Figure [21 shows the structure of the lower exciton en- 
ergy level as a function of the applied electric field Eqd 
in a CdSe nanocrystal 2 nm in radius. The lower exciton 
at zero electric field is 8-fold degenerate: 3 states with 
M = 1 and 5 states with M = 2. All these states are orig- 
inated from the ls-lS 3 / 2 electron-hole pairs. The electric 
field splits this level in two quartets: the lower one (I) be- 
longs to states with M z — ±1, ±2 while the higher level 
(II) in Fig. |21 corresponds to exciton wave function with 
M z = 0,0, ±1. The inset in Fig. |21 displays the matrix 
elements of the dipole operator, which determine the op- 
tical properties and will be discussed in the next section. 
TablelTTlillustrates the fraction contribution, \Cm(M z )\ 2 , 
of the dominant M-components in the expansion iJSJ to 
the excitonic levels I and II at Eqd = 500 kV/cm. It is 
worth to note that the states with M z = ±1 are derived 
from the M = 1 and 2 zero-field exciton wave functions, 
which are coupled by the electric field through remote 
states. The following group of exciton levels arises from 
the IS-IP3/2 electron-hole pairs, with a splitting pattern 
similar to that of the lower level. In general, all the 
M z = levels are doubly degenerate. It is important 
to remark that the above description is valid for larger 
nanocrystals in the strong confinement regime. 

The dependence of the binding energy upon Eqo , cor- 
responding to the I and II exciton levels are shown in 
Fig. |2| for several nanocrystal radii. 
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TABLE II: Fraction of the dominant M-component 
\C M (M Z )\ 2 to the expansion @ for the lower excitonic states 
I and II shown in Fig. [5] at the electric field intensity Eqd = 
500 kV/cm. 



Level 




I 


II 


E-E g 


390.7 meV 


393.4 meV 




\C M (M Z 


l^-fr action 


\C M 


(M z ) ^-fraction 


M\M Z 


(±2) 


(±1) 





(±1) 














0.01 


1 





0.74 


0.02 


0.96 0.25 
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0.99 


0.25 
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0.03 0.74 
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FIG. 3: Exciton binding energy versus electric field for the 
levels I and II shown in Fig. |5|for four QD radius (2,3,4, and 
5 nm). 




E QD (kV/cm) 

FIG. 4: Exciton energies with quantum number M z — vs 
electric field intensity for a CdSe nanocrystal 4 nm in radius. 

The energy levels with M z = 0, up to 100 meV above 
the lowest exciton, are plotted in Fig. 0] The states are 
labeled by the pure electron-hole pair contribution at zero 
field. Notably, no anti-crossing behavior is observed in 
levels ls-lS 3 / 2 and ls-lP 3 / 2 , and there is a weak anti- 
crossing between ls-lD 5 / 2 and ls-2S 3 / 2 in the region 100- 
150 kV/cm. The anticrossings should not be altered by 
the effect of the cubic terms in the hole Hamiltonian, 
although the absolute values will change. 



IV. FIELD INDUCED OPTICAL PROPERTIES 

In the optical experiments the intensity of the absorbed 
or emitted light is proportional to the exciton oscillator 



strength (squared absolute values of the dipole matrix 
element between the ground and excited states of the 
quantum dot). The field will modify the electron and 
hole wave functions which determine the exciton oscil- 
lator strength. As the field delocalizes the electron and 
holes in opposite directions, in principle the allowed ex- 
citonic transitions should decrease as the applied electric 
field increases. This is related to the quenching of the 
luminescence spectra^^ Moreover, the field breaks the 
inner spherical symmetry of the dot and forbidden exci- 
tonic transitions should appear in the optical spectra for 
large electric fields. 

The optical matrix element of the dipole operator for 
the exciton is given by 

D*, G = C , a(*)*D„, G> (11) 

a={n,N,l,L,f,F;M,Mz} 

where G represents the nanocrystal ground state and 
Y) a Q — i\e\h(a\p\G) / mo{E a — Eg) are the dipole ma- 
trix elements for the uncorrelated electron-hole pairs a 
(Eq. II U| l , which are given in Ref. |33j These dipole ma- 
trix element are different from zero only for M = 1 
and are proportional to e* M , where &m z =o = &z and 

ej|./ != ii = T (&x i i&y) /V% are the unit vectors in the 
spherical representation. This accounts for the optical 
selection rules in the dipole approximation: the only op- 
tically active electron-hole pairs have M z = 0,±1, and 
AI = 1 . According to the results of Table [I] the exci- 
ton level I in Fig. |2 is optically active for light polarized 
(with electric polarization vector e;) along e x or equiv- 
alently e/ || eM,=±i, while the state II is allowed for 
e; || e z || e M= =o or e/ || e x || e M ,=±i- The corresponding 
squares of the dipole matrix elements \D Z = e z ■ D^gI 2 
and \D X = e x ■ D^gI 2 for the states I and II are shown 
in the inset of Fig. [21 Note that at zero electric field the 
sum of | Ac | 2 over all M z — ±1 states equals to |-D Z | 2 , as 
required by the spherical symmetry of our Hamiltonian. 

Figure [S] shows the dependence of |-D Z | 2 upon the elec- 
tric field intensity for the lower excitonic states with 
quantum number M z = for CdSe nanocrystals of sev- 
eral radii. In the inset we present the same results 
shown in the figure but rescaled according to the laws 
(R /R)\D Z \ 2 and (R/R ) 3 E QD , where R = 2 nm. It 
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FIG. 5: Allowed dipole matrix elements, \D Z \ , of the lower 
excitonic states II as function of the applied electric field for 
CdSe nanocrystals of several radii. 



300 - 



' £ 200 

o 
Q 

Q~ 100 





R-4 nm 


IS- IS 3/2 


M=0 


_ls-2S 3/2 ls-lP 3/2 
}^UAD^ 


l s-l Fj/^T^^^^^^^ 







100 200 300 400 

E 0D (kV/cm) 

FIG. 6: Dipole matrix element \D Z \ 2 for a CdSe nanocrystal 
4 nm in radius as function of Eqd for different excitonic states 
with angular momentum projection M z = 0. 



the optical dipole is not quenched, while for large ra- 
dius the reduction in photoluminescence intensity should 
be significant. For example, for a field of 150 kV/cm the 
square dipole matrix element decreases approximately by 
66% for a QD 50 A in radius, while for a QD of 20 A, 
\D Z \ 2 is almost a constant for the range of the experi- 
mental values of the electric field that can be considered. 
The explanation of the above considered features lies in 
the interplay between the confined and electric energies, 
i.e. the kinetic energy and the external potential energy 
depend as R~ 2 and R, respectively. In Fig. the behav- 
ior of the oscillator strength \D Z \ 2 as a function of Eqd is 
shown for several excitonic states with angular momen- 
tum projection M z — 0. At Eqd = the exciton oscil- 
lator strength is diagonal and only transitions between 
electrons and holes with S- symmetry are allowed. The 
field breaks this selection rule and transitions with differ- 
ent symmetry are allowed, i.e. the electric field couples 
states with L, I ^ and the oscillator strength becomes 
different from zero for s — P, s — D, ... electron-hole 
pair transitions. The mixing effect, due to Eq/j, reduces 
the overlapping between electron and hole states with S- 
symmetry, allowing other electron-hole transitions. In 
Fig. El the transition IS-IS3/2 is the strongest one over 
the full range of electric field intensity. For higher tran- 
sitions, there are kinks at 300 and 480 kV/cm, which are 
due to avoided crossings between higher levels, e.g., be- 
tween IS-ID7/2 and IS-IF5/2 at 300 kV/cm (see Fig.0). 
Another interesting feature is the disappearance and the 
appearance of several transitions as the field is tuned. For 



25. 



3/2- 



shown 



example, the dipole of the transition Is 
as thick line in Fig. is strong at Eqd — and disap- 
pears at Eqd ~ 180 — 220 kV/cm. The opposite can be 
argued for the Is — I-D7/2 and Is — I-F7/2 transitions, for 
which the dipole matrix elements reach a maximum and 
at higher field intensities they have practical zero oscil- 
lator strength. Furthermore, the transition Is 
reappears at higher Eqd- 
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can be noticed that \D Z \ 2 scales almost as R, while the 
electric field scales as R~ 3 . The linear dependence of 
\D Z \ 2 upon R is a consequence of the Coulomb interac- 
tion, i.e. an excitonic effect. For sake of comparison in 
the Fig. [S] the free electron- hole calculations of \D Z \ 2 are 
shown by the triangles and circles for dot radii equal to 
2 nm and 4 nm, respectively. As can be seen in the fig- 
ure, the exciton effects are large, even for small radius 
of 2 nm. This means that the usual strong confinement 
approximation, where the Coulomb interaction is con- 
sidered as a small perturbation, breaks down for small 
CdSe nanocrystals. This effect can be explained due to 
the finite confinement barrier of the electron allowing the 
penetration of the exciton wave function in the surround- 
ing medium. Also, Fig. [S] indicates that for small radius 



Our calculation reproduces the magnitude of the Stark 
shift, which can be as large as 80 meV for inter- 
nal field intensity of 500 kV/cm in nanocrystals 4 nm 
in radius. This value agrees well in magnitude with 
the experimental observations of Ref. |8 if one consid- 
ers that the maximum of the optical gap indicates the 
zero of the internal electric field under an applied ex- 
ternal potential. The presence of an internal electric 
field and other features of single dot luminescence are 
currently attributed to trapped charges near the sur- 
face of the nanocrystal£2£i2LiSiSi Large dipole moments 
for the ground state and the LUMO in intrinsic CdSe 
nanocrystals were predicted by a set of pseudopotential 
calculations, 40 which can also account for the linear Stark 
shift in single dots. This intrinsic dipole moment is as- 
sociated with the lack of inversion symmetry in wurtzite 
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lattice and depends upon dot structural details and the 
dielectric response of the surrounding medium. The same 
calculation also reports null dipole moments for CdSe 
dots with zinc-blende structure, thus supporting our the- 
oretical model. 

Our results are qualitatively similar to those of Fu 41 
for InP QD, with the exception of non zero dipole mo- 
ment at zero field. Due to the inversion symmetry of 
our hole Hamiltonian (J2J, the hole eigenstates have defi- 
nite parity and zero dipole moment. A finite dipole mo- 
ment, caused by the lack of inversion symmetry in the 
zinc-blende structure, can be obtained if a linear term in 
k is included in the k • p Hamiltonian J2J. This effect 
is tiny in bulk semiconductors and is usually neglected. 
In nanocrystals the dipole may arise from mixing of the 
hole states IS3/2 and 1^3/2- The amount of mixing de- 
pends upon the ratio of the k linear term matrix element 
(oc 1/R) to the energy separation of the level (cx 1/R 2 ). 
Hence, the dipole should be small for small nanocrys- 
tals with zinc-blende lattice. For large nanocrystals the 
bulk regime is approached and the linear terms are again 
negligible. Although it is not possible to predict with cer- 
titude the dipole behavior for intermediate dot sizes, its 
absolute value should n ot g row. This is consistent with 
the calculations of Refs. l40l andl4ll 

An anomalous field dependence of emission lines of 
self- assembled quantum dots (SAQD) has been observed 
through micro-photoluminescence measurements^ where 
certain transitions lines of the luminescence spectrum dis- 
appear and reappear as the electric field is tuned. This 
feature could not be accounted for using a ID model of 
the quantum confinement £ The dipole matrix element of 
the state IS-2S3/2 (Fig. displays that behavior, which 
could be a property of the 3D confinement. However, a 
detailed calculation with the SAQD symmetry is needed 
to test this hypothesis. 

Let us consider two simple configurations for pho- 
toluminescence experiments. First, the emitted light, 
with wave vector n, is recorded along the Eqd direc- 
tion and its polarization vector e; is parallel to e x , (that 
is, e; || e x _L Eq/j). In this case the dipole element 
D x = e x ■ D^ g is not zero and, according to Table [IJ 
the states with M z = ±1 of the lowest excitonic level 
are optically active. Photoluminescence spectra should 
provide the Stark splitting presented in Fig. Second, 
the emitted light is observed perpendicular to the field 
(k _L Eqd). Here, we have two choices for the light po- 
larization, (i) e; || Eqc or (ii) e; _L Eq£>. In the case 
(i), e; || e z || etM z =o and the emitted light corresponds 
to the excitonic state II and the photoluminescence spec- 
trum presents only one peak. In the configuration (ii), 
both excitonic states I and II are activated and the Stark 
splitting of Fig. [3 appears. The case (ii) is illustrated 
in Fig. [7| for the absorption spectra of a single QD at 
different electric fields. In the figure, the quenching of 
the absorption lines and the split of the IS-IS3/2 exciton 
peak as the electric field increases are clearly observed. 

The above picture is slightly modified by the cubic 
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FIG. 7: Absorption spectra for a CdSe nanocrystal at dif- 
ferent electric field intensities. The light wave vector and its 
polarization e; are perpendicular to the electric field. 



anisotropy. At zero electric field, the 8-fold degener- 
ate IS-IS3/2 excitons belong to the representations 
(M = 1) and T>^ (M = 2) of group 0(3). It is remarkable 
that the Coulomb interaction does not remove the degen- 
eracy of T>y and T>% , and neither does it for higher levels. 
According to the compatibility tables for groups 0(3) and 
T d ^ = T 5 and V% = T 3 + T 4 . Hence, the M = 1 
triplet is not splitted. As only is dipole allowed, the 
selection rule M = 1 remains valid for the ground exci- 
ton state. Only exciton states ns-NS^/2 belong to 2)7/ 
and are dipole allowed in the spherical approximation. 
The rest of the states, except ns-NPi/ 2 , split in different 
levels that include T§ and should produce weak lines in 
the optical spectra. For non-zero electric field, the 0(3) 
symmetry is reduced to Coo„, and the irreducible repre- 
sentations correspond to the different values of \M Z \. The 
effect of the cubic terms depend upon the orientation of 
the crystal axes relative to the electric field. If Eqd is 
parallel to a 3-order axis, the symmetry is reduced to 
C% v . Hence, the excitons with M z = and M z = ±1 
belong to the Ti and T3 irreducible representations of 
C% v , respectively, and are not modified by the symmetry 
reduction. Excitons with higher \M Z \ can split and give 
rise to weak optical transitions. If the electric field is 
not oriented along a 3-order axis, all the degeneracies are 
removed and additional transitions should appear. How- 
ever, as the cubic anisotropy is induced fundamentally 
through the hole states, the splittings and intensities of 
the extra lines should be extremely small. 

The fine structure that we have revealed can be mod- 
ified by other effects present in real nanocrystals, such 
as shape asymmetry, crystal field (in wurtzite nanocrys- 
tals), and the electron- hole exchange interaction. 15 These 
effects lead to a splitting pattern of the 8-fold lowest ex- 
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citon into 5-levels, which is combined with the electric 
field-induced splitting. Non-adiabatic phonon-induced 
effects^ and dielectric mismatch-induced modifications 
of the Coulomb interaction^ contribute also to the mag- 
nitude of the splittings. For CdSe QDs with wurtzite 
structure the relative orientation of the crystallographic 
axes and the external electric field must have observable 
signatures in the optical spectra. These effects may be 
stronger than those produced by the cubic anisotropy. 
Since our approach ignores these effects our discussion is 
approximate. 

In summary, we have studied the influence of the elec- 
tric field on the electron and hole single particle states in 
CdSe nanocrystals, as well as on the exciton states and 
the optical properties. We have described the electric- 
field-induced quenching of the absorption and lumines- 
cence spectra and the importance of the exciton effects. 
We have found that the Coulomb interaction has a large 
influence on the strength of the optical transitions, even 
for small quantum dots. This fact is related to the 
penetration of electron wave function in the embedding 
medium, that partially breaks the strong confinement 
regime. Moreover, we have shown that for small QDs 
the dependence of the exciton binding energy upon the 
applied electric field is strongly reduced. For zero elec- 
tric field we have reported very general expressions for 
the solutions of the 4x4 hole Hamiltonian in spherical 
QDs. 
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APPENDIX A: MATRIX ELEMENTS OF THE 
STARK TERM 

The Stark term in the electron effective mass Hamilto- 
nian is a component of an irreducible spherical tensor of 
rank 1. The theorem of Wigner-Eckart and the reduction 
formulas for compound systems detailed in Ref. 01 allows 
us to write the matrix elements as 

(n'l'f'f' z \\e\E QD rcoad\nlff z ) 
= (-iy'-fW 2 +f\e\E QDy /(2f' + Wf + 1) 

x ( r i / \{ r i i \ (v i i \ 

\-f, fz J I 1 V 1/2 J I ) 
x A /(2Z' + l)(2/ + l) / R n ,i>(r)Rni(r)r 3 dr, (Al) 



where the terms inside {} and () are the Wigner's 6-j and 
3-j symbols, respectively. 

For the hole states we have found the expression 

(N'L'F'F' Z \-\e\E QD rcos6\ NLFF Z ) 
= (-l) F '- F *+V 2+F \e\E QD ^(2F> + 1)(2F+ 1) 

f F' 1 F \ ^ ( F' F 1 \ f K' 1 K\ 
X ^ -F' z F z ) 2_. \ K K> 3/2 / ^ J 

Xy/(2K' + 1)(2K+1) [ Rp K ,( r )R^ K (rydr.(A2) 



For the exciton states we found the expression 

(n'N'l'L'f'F 1 ; M' M' z \\e\E QD (r e cos6 e ~ r h cos6 h )\ nNlLfF; MM Z ) 

= \e\E QD (-l) M '-^S Mz .M,^ M M , I ™ z ) V(2M + 1)(2M' + l)L N , N >S L , L> S F , F ,(-^ F+M+f+f+1/2 

xV(2/' + l)(2/+l)(2/' + l)(2/ + l)| { {' i /2 } { M f f f] (o I) J Rn'v{r e )R nl {r e )rldr e 
-S n ^S l MS fJ/ (-l) M ' +f+ ^ 2+2F ^/(2F' + 1)(2F + 1) V(2tf' + l)(2tf+l) 

K,K' 

X { K K' 3/2 } { F F' / } ( ^ 0^)/ ^R'^J^&C^K^j- ( A3 ) 
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